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Abstract 

We consider both theoretically and experimentally self-organization process of quasi- 
equilibrium steady-state condensation of sputtered substance in accumulative ion- 
plasma devices. The self-organization effect is shown to be caused by self-consistent 
variations of the condensate temperature and the supersaturation of depositing 
atoms. On the basis of the phase-plane method, we find two different types of the 
self-organization process to be possible. Experimental data related to aluminum 
condensates are discussed to confirm self-organization nature of quasi-equilibrium 
steady-state condensation process. 
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1 Introduction 

Formation of various nanosystems in the course of the condensation process is 
known to be achieved due to widespread technologies such as the molecular- 
beam epitaxy, the metal-organic vapour-phase epitaxy, the liquid-phase epi- 
taxy, the electrolytic deposition etc. [1,2]. Characteristic peculiarities of these 
technologies are both quasi-equilibrium and steady-state to be as follows: 

• the condensation process occurs at sufficiently low supersaturation to pro- 
vide a proximity to phase equilibrium between condensate and depositing 
substance; 

• this supersaturation remains stable in the course of the time to guarantee 
steady-state conditions of the condensation process. 
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Due to above conditions adsorbed atoms can be coupled with growing con- 
densate surface only if the strongest chemical bonds are realized to minimize 
the free energy [3,4]. 

The proximity to the phase equilibrium is known to be governed by decrease 
of the supersaturation 



where n and n c are current and equilibrium concentrations of atoms above the 
growth surface. In the case of vapour-condensate systems, the latter of this 
concentrations is described with the empirical relation 



where Ed is the desorption energy, T is the growth surface temperature, ks 
is the Boltzmann constant; the temperature dependent parameter A(T) = 
exp(a + f3T + 7/T) is determined with a set of constants a, {3 and 7 being 
characteristics for a given substance [5]. 

In the case of volatile substances, the desorption energy is so small (Ed <C ksT) 
to expand large value of the equilibrium concentration (2) into a series. As a 
result, the supersaturation becomes less sensitive to temperature change that 
simplifies implementation of quasi-equilibrium steady-state condensation at 
relatively high deposition fluxes. 

For the vapour- condensate systems related to the molecular-beam epitaxy, 
low supersaturations are obtained only due to increased growth surface tem- 
perature at relatively weak deposition fluxes. That is why the molecular-beam 
epitaxy can be implemented mainly for those substances which have relatively 
high volatility. 

In chemically active medium- condensate systems, being the basis of metal- 
organic vapour-phase epitaxy, liquid-phase epitaxy and electrolytic deposition, 
the proximity to the phase equilibrium is stimulated additionally by reversible 
chemical processes. The latters induce reverse transitions of adsorbed atoms 
which are weakly coupled with growth surface into environment. In Eqs. (1) 
and (2) this is expressed in decreasing value E d of the desorption energy down 
to an effective value that increases the equilibrium concentration n c stimulat- 
ing proximity to equilibrium. From this point of view, above chemical meth- 
ods make it possible to reach low supersaturations at condensation of weakly 
volatile substances that gains its advantage over the molecular-beam epitaxy. 

The present work is devoted to consideration of the system low temperature 
plasma- condensate which enables to reach a proximity to phase equilibrium by 
means of heating growth surface with the help of plasma stream, both momen- 
tum and energy transfer from plasma particles to adatoms, and their partial 
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thermal accommodation. Due to the latter two mechanisms the desorption 
energy E d gets an effective value [3,6,7] 



E = E d -5E 



(3) 



with a stochastic addition 5E characterized by the mean value E and the dis- 



persion a 2 = y5E — Ej (the overline means averaging over plasma particle 
collisions). Thus, in analogy with the chemical systems the low-temperature 
plasma-condensate systems make it possible to stimulate proximity to equi- 
librium at condensation of weakly volatile substances. However, making use 
of the plasma arrives at the problem to fix steady-state conditions during the 
condensation process. 

For all mentioned technologies, there arises a common range of problems to fix 
steady-state regime of the condensation process at vanishing supersaturations. 
The main point is that such crucial technological parameters as the growth 
surface temperature T (hence, the equilibrium concentration n e = n e (T)) and 
the depositing flux J are mutually independent and are regulated by different 
power sources. Also, under quasi-equilibrium conditions the condensing flux J c 
is obeyed the inequality J c <C J. Hence, even slight relative variations A J/ J of 
the depositing flux can arrive at considerable changes in condensation kinetics 
at condition A J ~ J c . (The same effect can cause variation of the equilibrium 
concentration n e (T) at slight fluctuations of the growth surface temperature 
T.) In the absence of self-organization, these problems are solved by artificially 
created feedback between condensation kinetics and depositing flux by means 
of control systems that arrives at considerable rise in the cost of technologies. 

Therefore, the aim of this Letter is two-fold: i) elaborate a technological de- 
vice, whose operation within self-organization regime ensures extremely low 
steady-state supersaturations; and ii) build up a theoretical model to describe 
above self-organized process. In correspondence with these aims, the outline of 
the paper is as follows. In Section 2, we state physical backgrounds of opera- 
tion of the accumulative ion-plasma devices. On the basis of these conceptions, 
we consider theoretically the self-organization of extremely low stationary su- 
persaturations in Section 3. Experimental validation of the self-organization 
scheme is given in Section 4, and Section 5 concludes our consideration. 



2 Physical backgrounds of operation of the accumulative ion-plasma 



Principle components of our sputtering device are the magnetron sputterer 
in combination with the hollow cathode which operate under high pressure of 
highly refined argon (P Ar = 10-=-30Pa) (see Fig.l). Such a construction allows 




device 
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Fig. 1. (Colour online) Scheme of axisymmetric ion-plasma device (1 - anode, 2 - 
thermalization volume of sputtered atoms, 3 - erosion zone, 4 - magnet system, 5 - 
hollow cathode, 6 - region of circular mass transfer, 7 - layer of adsorbed atoms, 8 
- condensate, 9 - substrate, 10 - chiller); plot on the left depicts both temperature 
and concentration distributions near growth surface 7. 

to stabilize and increase the discharge current due to both hollow cathode 
and the magnetron effects [8,9], the last of which localizes the erosion zone 
in addition. It is worthwhile to stress that almost all atoms of inert gas and 
sputtered substance become ionized when get inside the hollow cathode [8,9]. 
Under increased pressure, plasma particles collide so frequently and intensively 
to average their energy effectively. These collisions decrease the dispersion a 
of the effective desorption energy (3) that, in turn, stabilizes condensation 
process. 

The self-organization process starts with accumulation of a sputtered sub- 
stance near the growth surface. According to Ref. [10], atoms sputtered by 
magnetron lose their energy (in other word, thermalize) due to collisions of 
particles of high-pressure plasma, so that their further movement becomes dif- 
fusive. Related thermalization length is estimated as I oc {Tat + AT) /TatPat 
with Tat and AT being the temperature of unheated working gas and its 
growth under discharge. Hence, the thermalization length is reduced with in- 
creasing the pressure Pa t and decreasing discharge power that reduces the tem- 
perature growth AT. As show estimations [9,10], the flux of non-thermalized 
atoms in direction of the anode can be considered as negligible at the pres- 



4 



sure Pat = 20 -j- 40 Pa and the discharge power less than 100 W, when the 
thermalization length is subject to the condition I < L/4 with respect to the 
target-anode distance L depicted in Fig.l. 

As it is seen from this figure, the flux J s of sputtered substance condenses 
mainly on the surface adjacent to the erosion zone with the flux J q and diffuses 
partially into the hollow cathode with the flux J d . (It should be noted that at 
the initial time of the device operation, there is only inert gas plasma inside 
the hollow cathode.) With power up, the concentration of sputtered atoms n 
at the hollow cathode entry is determined by the fluxes J e and Jp, the first 
one being drift flux due to the presence of electrical field near the inlet, and 
the second one being caused by pressure difference between interior and entry 
of the hollow cathode. As a result, variation of the atom concentration no at 
the hollow cathode inlet is governed by the equation 

h oc (J s -J d - J q ) - (J e - J P ) (4) 

where and hereinafter the point above a symbol denotes the time derivative. 
As the hollow cathode volume is rather small, the drift components J e and 
Jp become equal quickly and then sputtered substance penetrates into the 
hollow cathode due to diffusion only. Because both components J s and J q are 
much more than the diffusive flux J d , one can suppose its variation affects 
slightly on the concentration n . At the steady-state condition J d = J s — J q , 
diffusion process determines accumulation of thermalized atoms in the interior 
of the hollow cathode to increase their concentration over a thermodynamic 
threshold n e needed to start condensation of the sputtered substance on the 
substrate P~l 

To understand the lowering of steady-state supersaturation in the course of 
the condensation process, let us assume the growing surface has reached such a 
temperature, when inequality n > n c starts to be satisfied due to accumulation 
of sputtered substance. Then, the supersaturation (1) takes a positive value 
to initiate the barrier nucleation process that increases the growth surface 
temperature due to the atom thermal accommodation. On the other hand, 
the flux of depositing atoms decreases due to rupture of circular fluxes that, 
in turn, reduces the supersaturation to minimal values providing condensation 
process. In the following section, we shall show that self-consistent variations 
of the surface temperature and the supersaturation is a cornerstone of self- 
organization process ensuring extremely low steady-state supersaturations. 



It is worthwhile to note, under plasma influence, main part of deposited atoms 
are reevaporated to be next ionized again and returned onto this surface under elec- 
tric field [8] . Such a circular mechanism of mass transfer accumulates additionally 
depositing atoms in immediate vicinity of the growth surface. 
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3 Theoretical consideration of self-organization of steady-state quasi- 
equilibrium condensation 



Let us start with consideration of main reasons for variation of the surface 
temperature in the course of the ion-plasma device operation. First, each of 
plasma ions transfers to the surface the average energy k B (T 2 — T) deter- 
mined by the temperature difference of plasma ions, T 2 , and adsorbed atoms, 
Ti, respectively (see the temperature distribution depicted in the left panel of 
Fig.l). Then, the energy delivered per unit time from the plasma to the unit 
surface can be written in the form 

E 1 = X 0(T 2 -T). (5) 

Here, x is a parameter being the production of the Boltzmann constant and 
the plasma flux falling down the growth surface, 9 = (T 2 — Ti)/(T 2 — T) ~ 
1 — Ti/T 2 = const is thermal accommodation coefficient of adsorbed atoms 
rewritten with accounting the condition T T 2 . Second, the rate of energy 
transfer to the growth surface due to thermal accommodation of condensed 
atoms 

E 2 = k B (T 2 -T)J c (6) 
is determined by the temperature difference of plasma and condensate accom- 
panied with the condensing flux 

J c =-(n-n e ). (7) 
r 

Here, S and r are characteristic length and time of the circle motion of con- 
densing atoms near the growth surface. Finally, the energy removal to the 
chiller per unit time by means of thermal conductivity is determined by the 
obvious expression 

£ 3 ~^(T -T) (8) 

where To stands for the chiller temperature, d and 77 are the total thickness 
and effective thermal conductivity of two-layer system consisting of condensate 
and substrate. 



Combination of the equations (5) - (8) with obvious definition of the variation 
rate of the surface temperature 



cT = E 1 + E 2 + £ 3 
arrives at the differential equation 



(9) 



cT = 



\0T> ■ + 
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T 



T 2 £ 



(10) 
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with c being the surface heat capacity of the growth surface. The equation (10) 
shows the time dependence of the surface temperature is determined by the 
following set of values: plasma parameters Xi & an d t; condensate-substrate 
parameters 77, c and d; temperatures characterising condensate (T), plasma 
(Ti,T 2 , 9 ~ 1 — Ti/T 2 ) and chiller (To); equilibrium concentration n e defined 
with Eq.(2); and finally, supersaturation £ given by Eq.(l). Remarkably, both 
square brackets in Eq.(10) depend on the temperature T only through the 
equilibrium concentration (2), whereas the last term in right-hand side con- 
tains the factor T itself. Moreover, the last terms in both square brackets are 
proportional to the supersaturation £, whose time dependence must be found 
next. 

At finding the equation for the concentration n, one should take into account 
that condensation occurs onto the inner surface of the hollow cathode with 
large area S, whereas the diffusive flux penetrates toward the cathode interior 
throughout the inlet hole with small area s. Then, the rate of the concentra- 
tion n variation is determined by the difference between incoming diffusive 
and outgoing condensing fluxes so that, with accounting that accumulation 
of substance takes place inside the region with thickness 5 near the growth 
surface, one obtains 

» = ] (| * - j.) . (id 

According to Eq.(ll) the steady-state regime is reached under the condition 

sJ d = SJ C (12) 

which arrives at the relation J c < J d due to the non-equality s < S. The 
condensing flux is determined by Eq.(7), as the diffusive component is given 
by the Onsager-type expression 

Jd~-y(n-n ) (13) 

with effective length A, along which the concentration varies from n to n, and 
the temperature dependent diffusion coefficient D ~ T^ 2 [10]. Taking into 
account Eqs. (2), (7), (11) and (13), we yield the final form of the equation of 
motion for the supersaturation (1): 

f +B <™ 1+ f>*=SiM; + S£)< (14) 

where B{T) = (3 — ^ + ( Ed ^f)--y an j a ma g n itude £ = ( n o — n c) /n e relates 
to the inlet hole of the hollow cathode. Similarly to the differential equation 
(10), the first term in right-hand side of Eq. (14) depends on the temperature 
only through the equilibrium concentration (2), whereas the second term is 
proportional to the supersaturation with a constant factor inside the paren- 
theses. 
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The system of the ordinary differential equations (10) and (14) describes self- 
consistent variations of the surface temperature T and the supersaturation £ 
during the self-organization process of extremely low steady-state supersatura- 
tions. The form of possible solutions of these equations is seen from the phase 
portraits depicted in Fig. 2 where the set of parameters is selected applicable 
to the aluminum. Analysis of the phase portraits arrives at the conclusion 
that regime of the self-organization is governed specifically by the ratio s/S. 
In the case of high values s/S (Fig.2a), with the temperature growth the rep- 
resentative phase path 1 goes through high supersaturations, whereas in the 
opposite case (Fig. 2b) negative values £ precede to the steady state related to 
the node. From the physical point of view, such a behaviour is explained by 
that accumulation of the sputtered substance is in advance of heating of the 
growth surface, in the first case, whereas at small ratios s/S slight accumula- 
tion rate enables positive supersaturations £ only near the steady state. From 
the technological point of view, the second regime is much more preferable in 
case of little deposition time. 

Under the steady-state conditions T = £ = 0, the system of Eqs. (10) and 
(14) shows that the stationary supersaturation increases monotonically with 
growing effective desorption energy E, so that extremely low supersaturations 
£ <C 1 are achieved only when E < 0.1 eV. Usually, the bare desorption energy 
for metals is estimated as E^ ~ 0.4 eV [5], and the proximity to the phase 
equilibrium is provided by the plasma influence that reduces the bare value 
Ed to the effective energy E in accordance with Eq.(3). 



4 Experimental validation of self-organization 

We have considered above the self-organization process of quasi-equilibrium 
steady-state condensation of sputtered substance in accumulative ion-plasma 
devices. At the same time, both numerical simulations [11,12,13,14] and exper- 
imental investigations [3,14,15,16,17,18] display a secondary self-organization 
process. In the course of such a process, three-dimensional micro- and nanos- 
tructures are arranged on the substrate to possess complicated architecture 
with narrow distributions over sizes and forms of structural elements. In con- 
trast to the self-organization of the system plasma-condensate, the secondary 
self-organization involves in the regime of the self- organized criticality driven 
by fluctuations [19]. Really, the result of the primary process is that two phase 
system plasma-condensate is staying near stationary equilibrium, but strongly 
above a critical point, so that a phase transition is thermodynamically impossi- 
ble. However, appearance of active centers, playing the role of space distributed 
fluctuations, initiate the condensation process according to the following sce- 
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Fig. 2. (Colour online) Phase portraits presenting solutions of Eqs. (10) and (14) 
at s/S = 1/4 (a) and s/S = 1/10 (b) (the rest of parameters is as follows: 
D = 10 5 cm 2 s-\ E = 2.674 • 10~ 20 J, d = 0.1cm, A = 1cm, 5 = 0.1cm, r = 10~ 6 s, 
n = 10 8 cm- 3 , 6 = 0.95, x = 7.4 ■ 10- 3 J- 1 ci7»- 2 a -1 J T 2 = 6000K, T = 300-FsT, 
c = 3 • 10~ 3 J~ , r\ = 6 • lO -3 !^ cm. -1 ^ -1 (for glass). The node corresponds to 
T = 898K, £ = 1.1 • 10" 4 , J c = 1.1 • lO^m"^- 1 (a) and T = 898K, f = 7.4 • 10~ 5 , 
J c = 7.4 • K^cm -2 * -1 (b). 

nario]" 2 ! 

2 The condensation mechanism is beyond the scope of the equilibrium thermody- 
namics to be defined by both kinetics [4] and synergetics [19]. 
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In the course of the deposition process, adsorbed substance is embedded atom- 
by-atom into the growing surface on active centers with the highest energies of 
chemical bonds. Among such centers, one can point out the mono-step bends 
on the growing crystal surface, inhomogeneities on the atomic-rough surface, 
regions of crystallite joining and so on [5]. Due to the reduced surface density 
of the active centers, the chemical bonds spectrum should be appreciably dis- 
crete. That allows for one to separate those centers whose energies of chemical 
bonds Ei, i — 1, 2, . . . is lower than a critical threshold E c fixed by external 
conditions]^ On these centers, the atom-by-atom formation of the condensate 
becomes hardly probable because of subcritical energies of chemical bonds 
Ei < E c . By definition, the threshold E c relates to the desorption energy Ed 
in the equilibrium concentration (2) at zero supersaturation (1) and given con- 
centration n and temperature T. According to Eqs. (1) and (2), the critical 
energy of chemical bonds 

(15) 

increases with tending to the phase equilibrium due to concentration decrease 
and temperature growth. 

In our experiment, we investigate the structure of the aluminium conden- 
sates, produced in proximity to equilibrium in the plasma-condensate system. 
The scanning electron microscopy images of the ex-situ grown condensates 
are shown in Fig. 3. These condensates were obtained on glass substrates by 
means of the accumulative ion-plasma system described in Section 2 at the 
discharge power 1.8W during the deposition time t = 9h and at both pressures 
20Pa and 15Pa of the highly refined argon. It is the argon pressure as tech- 
nological parameter, that controls supersaturation, or degree of proximity to 
equilibrium, at stable discharge power. Thus, if argon pressure decreases from 
20 to 15 Pa, plasma particle concentration decreases as logical result, and, 
therefore, their average energy increases. This reduces the effective desorption 
energy and accordingly to Eqs. (1) and (2) reduces supersaturation as well. 

According to Fig. 3a, at increased pressure 20Pa, there are condensed weakly 
bound crystals whose habitus is fixed by crystallography planes akin to (210) 
because of the the conditions £(531) < E c < E( 2 io) inherent in energies E^m 
of atoms in half-crystal position on the (/i/c/)-planes of fcc-metals [20]. It is 
worthwhile to stress the structure depicted in Fig. 3a can be statistically ho- 
mogeneous up to tens of micrometers in thickness. Such a picture of the con- 
densate formation is caused with the sequence of repeated nucleation on the 
active centers of the growth surface being usually regions of crystallite joining. 

With reduction of the working gas pressure down to 15Pa, at remaining sta- 

3 Energies of chemical bonds Ei are considered as their absolute values. 



E r = k B T In 



nk B T 
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Fig. 3. (Colour online) Scanning electron microscopy images of ex-situ grown alu- 
minium condensates at layer-by-layer (a) and normal (b) crystal growth (the arrows 
on the lower panel point out the single crystals retaining faceting of the previous 
structure). There are used the glass substrates at the argon pressures 20Pa (a) and 
15Pa (b), the discharge power the deposition time 9h, the area ratio S/s ~ 4 

(the condensate thickness is of order 12 4- 14fj,m). 
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ble of the rest of technological parameters, the supersaturation decreases 
approaching the phase equilibrium of the system plasma-condensate. This 
means increasing the critical energy of chemical bonds E c which becomes 
higher the energy £'(210) of atoms on the plane (210), but lower the energy 
E r of chemical bonds on atomically rough surface. As a result, the conditions 
-E'(2io) < E c < E r are fulfilled to change the mechanism of atom attaching to 
growth surface. Above assumption explains the transition from tangential to 
normal growth of round shaped crystals with atomically rough surface shown 
in Fig. 3. Proximity of normal and tangential mechanisms of the crystal growth 
is confirmed by the fact that in Fig. 3b several crystals labeled with arrows have 
faceting appropriate to tangential growth along the (210) plane. It should be 
stressed as well that, due to quasi-equilibrium conditions of the condensation 
process, there is no secondary nucleation in the regions of crystallite joining 
what governs formation of the columnar structures [21]. 

Thus, the comparison of the electron microscopy images depicted in Figs. 3a 
and 3b allows for one to conclude that the microstructures related are obtained 
as result of different mechanisms of crystal growth. Obviously, above transition 
from layer-by-layer to normal crystal growth may be possible at two-fold condi- 
tions: first, the system plasma-condensate should be extremely near the phase 
equilibrium (to increase the critical energy E c ); second, condensation process 
should be strongly steady-state (to fix a location of the threshold E c in the 
chemical bond spectrum). As is shown in Section 3, both conditions pointed 
out are provided due to the self-organization of quasi-equilibrium steady-state 
condensation within device presented in Section 2. 



5 Conclusion 

In contrast to widespread technologies, we propose here original method to 
realize quasi-equilibrium steady-state conditions of condensation. Character 
peculiarity of this method is that condensate surface grows atom-by-atom on 
active centers having the strongest chemical bonds. We show such regime is 
ensured by the accumulative ion-plasma device whose advantage in comparison 
with already existing systems consists in that it does not require cumbersome 
control system. The self-organized regime of the device operation is confirmed 
both theoretically and experimentally. 

Above regime is shown to be achieved by means of two principle components 
of the device to be the magnetron sputterer and the hollow cathode. Due 
to diffusion, accumulation of condensed atoms occurs in the interior of the 
hollow cathode up to quasi-equilibrium concentration whose value remains 
constant during condensation process. Heating growth surface with the help 
of both plasma stream and decreasing the effective desorption energy causes 
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extremely low supersaturation near condensate overall. 

The physical reason to reach the quasi-equilibrium steady-state conditions 
is that condensation process evolves under scenario of the self-organization. 
On the basis of the phase-plane method, it is shown the self-organization 
is governed by self-consistent variations of the surface temperature and the 
supersaturation. 

To confirm that quasi-equilibrium steady-state condensation process is evolved 
under scenario of the self-organization, we investigate the structure of the alu- 
minium condensates obtained at different technological conditions. Related 
scanning electron microscopy images show the microstructures obtained as re- 
sult of the transforming layer-by-layer crystal growth into normal one. That, 
in turn, may be possible if the system plasma-condensate is near the phase 
equilibrium and condensation process is strongly steady-state. Just that very 
conditions are provided due to the self-organization of low steady-state super- 
saturations within device proposed. 
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